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Light hadronic physics using domain wall fermions in quenched 
lattice QCD* 

Matthew Wingate'^ 

'^RIKEN BNL Research Center, Brookhaven National Laboratory, Upton, NY 11973, USA 

In the past year domain wall fermion simulations have moved from exploratory stages to the point where 
systematic eflects can be studied with different gauge couplings, volumes, and lengths in the fifth dimension. 
Results are presented here for the chiral condensate, the light hadron spectrum, and the strange quark mass. 
We focus especially on the pseudoscalar meson mass and show that, in small volume, the correlators used to 
compute it can be contaminated to different degrees by topological zero modes. In large volume a nonlinear 
extrapolation to the chiral limit, e.g. as expected from quenched chiral perturbation theory, is needed in order to 
have a consistent picture of low energy chiral symmetry breaking effects. 
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1. INTRODUCTION 

The RBC Collaboration has recently reported 
results based on quenched QCD simulations with 
domain wall fermions on lattices of various vol- 
umes and spacings IT] . This talk summarizes our 
understanding of the behavior of the pseudoscalar 
meson mass as a function of input quark mass, 
scaling of the nucleon-rho mass ratio, and strange 
quark mass. Another set of quenched simulation 
results with domain wall fermions have been re- 
ported by the CP-PACS Collaboration §. 

2. CHIRAL CONDENSATE 

A Dirac operator D which has a well-defined 
index should lead to a chiral condensate which 
satisfies the following Banks-Casher relation (for 
finite volume) 



-{qq) 



where A^ are eigenvalues of the massless Dirac op- 
erator and v is the index of D. Note that these 
zero modes are chiral and correspond to units of 
topological charge in the continuum. In full QCD 
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Figure 1. Chiral condensate vs. mass for volumes 
8^ X 32 (circles) and 16^ x 32 (squares) at /3 = 5.7, 
L, = 32. 



simulations configurations which support eigen- 
values equal to zero are suppressed during Monte 
Carlo updating due to the fermion determinant in 
the Boltzmann weight. This suppression is absent 
in the quenched approximation. Consequently, a 
prominent 1/m divergence should appear in (qq) 
for quenched simulations in finite volumes. Only 
when one uses a fermion discretization which ad- 
mits a non-zero index can one observe this diver- 
gence, so the first signal was seen with domain 
wall fermions y,^ . 

Fig. ffl shows (qq) vs. the input quark mass m/ 
on two volumes at /3 = 5.7 (plaquette gauge ac- 
tion). As expected from Eq. |l| the divergence 
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Figure 2. Monte Carlo evolution of correlators 
on a 8^ X 32 lattice at /? = 5.7 with Ls = 48 and 
m/ = 0.0. 



is much more severe for the smaller volume, ap- 
proximately (1.6 fm)^, than for the larger volume, 
Ri (3.2 fm)"^. The lines are fits to the form 

- (ot) = —; h ao + ainif . (2) 

(qq) has a similar divergence on a (1.6 fm)'^ box 
at /3 = 6.0 §. 

3. PSEUDOSCALAR MESON 

Given that the chiral condensate in finite vol- 
ume is sensitive to topological zero modes when 
using domain wall fermions, we turn to corre- 
lators which couple to the pseudoscalar meson, 
viz {A^{x)A^{0)) and (7r''(2;)7r''(0)}. Looking at 
the spectral decomposition of these correlators |1|] 
leads one to expect the leading zero mode effects 
to be {A^{x)A^{0)) oc 1/m and (7r"(x)7r'^(0)) oc 
1/m? . Indeed time histories of these correlators 
(Fig. g) show the latter correlator to be much 
more singular than the former. Also plotted in 
Fig. is the connected scalar-scalar correlator 



Figure 3. Fits to pseudoscalar meson mass 
from (7r'^(a;)7r'^(0)) (D), (Ag(a;)Ag(0)) (o) and 
(7r'^(a;)7r''(0)} + {a{x)a{{)))c (O) on a 16^ x 32 lat- 
tice at /3 = 6.0 with Ls — 16. The solid line is 
a linear fit to the circles with 0.01 < ruf < 0.04. 
The asterisk marks m,-cs and the inset magnifies 
the data at small mf. 



which has the same leading zero mode effects as 
(7r'^(a;)7r'^(0)), both according to the spectral de- 
composition and as seen in the data. 

The singular behavior of {tt'^ {x)tt°- (0)) for small 
mass leads to a contamination in the extrac- 
tion of the meson mass. The mass computed 
from this correlator is significantly larger at small 
nif than the mass computed from {A^{x)Aq{0)), 
both at /3 = 6.0 (Fig. |) and at /3 = 5.7 (Fig. 
^). An attempt can be made to subtract the 
zero mode effects by extracting a mass from 
(7r''(a;)7r''(0)) + {a{x)a{0))c. At small m/ the fact 
that the mass from {Aq{x)Aq{0)) is larger than 
from the "subtracted" correlator could be due to 
1/m zero mode contributions in {Aq{x)Aq{0)). 
The subtracted correlator receives increasing con- 
tribution from the isovector scalar meson as m/ 
increases, explaining the deviation of the dia- 
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Figure 4. Fits to pseudoscalar meson mass from 
correlators denoted as in Fig. gon a S'^ x 32 lattice 
a.t f3 — 5.7 with Lg — 48. The solid line is a linear 
fit to the circles with 0.02 < m/ < O.f . 



monds for larger values of irif in Fig. y. 

The zero mode effects disappear within the sta- 
tistical resolution when the volume is increased. 
Fig. ^ shows the pseudoscalar meson mass ex- 
tracted from all three propagators on a lattice 
with (3.2 fm)^ spatial volume at /3 = 5.7. The as- 
terisk in that figure denotes the residual mass, 
rriics — 0.0072(9), computed from the ratio of 
the mid-point pseudoscalar density to the sur- 
face pseudoscalar density WH . Although it is not 
clear on the scale of the plot, the linear extrapola- 
tion of (am^)^ to zero in mf misses ?Tiros by over 
2(7. Quenched chiral perturbation theory suggests 
that for vanishing quark mass m, rn^ oc 771^''^+'' 
IgI. In order to simplify the fit in the region of 
mf where the data points are, one uses 

(arTiTr)^ ~ aa[mf + ai)(l — (51n(r7i/ + ai)) . (3) 

The sohd hue in Fig. | gives ai = 0.0073(10), 
in good agreement with m-^cs, and 5 = 0.07(4). 
Although the x^/dof (with errors estimated from 



Figure 5. Pseudoscalar meson mass squared vs. 
m/ on a Idi^ x 32 lattice sX (3 — 5.7 with Lg = 48. 
Symbols used as in Fig. ^ and the asterisk marks 
m-ics- Solid line is a quenched chiral log fit to 
circles with 0.0 < mf < 0.08, and the dashed line 
is a linear fit to the same data. 



jackknifing) decreases from 4.3±2.6to3.6±2.4 
between the linear and logarithmic fits, it is not 
compelling enough to conclude firmly that the lat- 
ter is preferable. On the other hand, our theoret- 
ical prejudice is for a fit which extrapolates to 
(am-^)'^ = at ruf = —m-^cs- 

The important conclusion from this study is 
that topological zero modes can contaminate 
quantities used to compute observables in small 
volume simulations. This interesting effect is a 
consequence of the index of the domain wall Dirac 
operator and, and it vanishes for increasing vol- 
ume. 



4. SCALING OF m^/mp 

Vector meson and nucleon masses have been 
computed as functions of nif on lattices sX (3 ~ 
5.7, 5.85, and 6.0 on approximately equal spa- 



tial volumes. Fig. ga shows ruN /rnp vs. a^ where 
the nucleon and rho masses have been extrapo- 
lated linearly to the chiral point ruf = —m-ccs- 
Although the scatter of the data points is only 
at the 1(7 level, it is hard to draw a firm conclu- 
sion regarding their scaling behavior. The fit to 
Co -I- 020? has a x^/dof of 2.4. 

Given that the data for the pseudoscalar data 
hint at nonlinearities in the chiral extrapolations, 
one might wonder what role these would play in 
Fig. ||a. The smallest rrif in the (3 = 5.7, 5.85, 
and 6.0 simulations correspond to mT^/mp = 0.38, 
0.57, and 0.46, respectively. Therefore the differ- 
ent simulations would have different sensitivities 
to any deviation from linearity in the chiral limit. 
One way to test scaling while avoiding chiral ex- 
trapolation is to compare rriAr/mp at a value of 
m-^/mp where simulations have been done. In 
Fig. @b we plot rriN /rup interpolating the data to 
mT^jrap = 0.61. These data fit cq 
xVdof of 0.88. 
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5. STRANGE QUARK MASS 

The calculation of the strange quark mass using 
this data was presented at DPF 2000 0. Using 
masses and matrix elements from |1| the strange 
quark mass is computed defined through both 
the vector and axial Ward-Takahashi identities 
(VWTI/AWTI). The renormahzation factors are 
computed in the RI/MOM scheme [|j|]. Wc find 
at /3 = 6.0 



m 



MS 



(2 GeV) 



110(2) (22) MeV VWTI 
105(6)(21) MeV AWTI 



(4) 



using the K mass to fix mf for the strange sector. 
The first error is statistical, and the second is 
the systematic error due to using perturbation 
theory at 2 GeV to match to the MS scheme. 
At /3 = 5.85 we quote mp(2 GeV) = 100(5) (20) 
MeV from the AWTI. Using the K* instead to set 
m^f increases mf^{2 GeV) by w 20% at /3 = 6.0 
and by « 35% at /3 = 5.85. 
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Figure 6. Nucleon-rho mass ratio on spatial 
volumes « (1.6 fm)'^. Squares denote domain wall 
results at /3 = 6.0, 5.85, and 5.7. The chiral limit 
ruf + TOros —^ has been taken in (a), and the 
asterisk marks the real world value, (b) shows 
mjsf /rup at mT^/mp — 0.61. Dashed lines show fits 
to Co -I- 020? and fancy squares the extrapolated 
a = value. 
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